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The recently proposed scheme for direct sampling of the quantum phase space by 
photon counting is discussed within the Wigncr function formalism. 



First complete experimental characterisation oLthe quantum state of a single 
light mode was demonstrated by Smithey et auH In their work, numerical to- 
mographic algorithms were applied to reconstruct the Wigner function from 
homodyne statistics. Recently, a novel scheme for measuring the quantum 
state of light by photon counting has been proposedB In contrast to tomo- 
graphic techniques, the newly proposed method allows one to measure qua- 
sidistribution functions directly at a selected point of the phase space. In this 
contribution we show that the Wigner phase space representation provides a 
clear and intuitive interpretation of this scheme. 

The principle of direct sampling of the quantum phase space by photon 
counting is to measure the photon statistics {p n } of the signal field as, su- 
perposed on a probe field ap by means of a beam splitter with the power 
transmission T. The measured photocount statistics is used to calculate the 
alternating series X)^Lo( — ty n Pni which is given by the expectation value of 
the following normally ordered operator expressed in terms of the signal and 
the probe fields: 

ft = : exp[-2(v / Ta| - Vl - Ta P )(Vfa s - Vl - Top)] : . (1) 

We will interpret this observable using the quantum mechanical phase 
space formalism introduced by Eugene Wigner. For this purpose we will rep- 
resent II as an integral of the product of two quantities depending only on the 
signal or the probe mode operators: 

n = lj d2 P '■ exp[-2(Vr/r - ap)(Vrp - a P )} ■. 

x : exp[-2(y/T=T0* - a f s )(VT^T(3 - a s )] : . (2) 

As the P mode is to be used as a probe for the S mode, we assume that both 
the modes are not correlated. Evaluating the quantum expectation value of 
the above integral representation yields the following formula: 



<n) = T) I d 2 /3 W s (0)Wp(y/T/(l-T)/3), (3) 
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Figure 1 : Phase space interpre- 
tation of the observable A mea- 
sured in the photon counting 
experiment. 



which shows that the measured quantity is simply the phase space integral of 
the product of the signal and the probe Wigner functions Ws{fi) and Wp((3). 
Thus, the quasidistribution function in the signal phase space is "sampled" by 
the probe Wigner function. Of course, the area of the "patch" of the probe 
Wigner function must be at least that imposed by the uncertainty principle. 
The essential advantage of the newly proposed scheme is that the param- 
eterisation of the probe Wigner function in Eq. (3) is rescalcd by the factor 
yjT/(l - T), which may take any positive value for < T < 1. When T > 1/2, 
the probe Wigner function becomes effectively "squeezed" in all directions, and 
consequently the resolution of the phase space measurement goes beyond the 
uncertainty relation limit. This feature is depicted schematically in Fig. 1. 
Let us stress that the rescaled Wigner function does not appear physically in 
the setup; it is only a tool for the phase space interpretation of the measured 
observable. 

An important class of probe states which can be easily generated in a 
laboratory are coherent states \a). In this particular case the outcome of the 
measurement can be expressed using an s = — (1 — T)/T ordered quasidistri- 
bution function of the signal mode: 



In the limit T — > 1 the ordering parameter approaches zero, which corre- 
sponds to the direct measurement of the Wigner function of the signal field. 
The Wigner function is determined at a point defined by the amplitude and 
the phase of the probe coherent state. By changing these two parameters, 
the complete Wigner function can be scanned point-by-point, without using 
complex numerical algorithms. 
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Figure 2: The one— photon Fock state quasidistribution function reconstructed from 1000 
Monte Carlo events at each phase space point (a) without and (b) with compensation for 
the non— unit detector efficiency r] = 80%, in the limit T — * 1. 



This result is easily understood within the Wigner function formalism. In 
the limit T — > 1 the rescaled Wigner function of the coherent probe approaches 
the shape of the delta function, and therefore the integral in Eq. (3) picks up the 
value of the signal Wigner function at a single phase space point. Additionally, 
the rescaling moves the probe Wigner function towards the phase space origin, 
which is reflected by the decreasing factor multiplying a in Eq. (4). Therefore, 
large intensity probe fields have to be used to scan the required region of the 
phase space. 

In a realistic setup, there is a limitation on the resolution of the phase space 
sampling imposed by the non-unit quantum efficiency r\ of the photodetector. 
Furthermore, rigorous statistical analysis shows that an attempt of numerical 
compensation for this deleterious effect fails due to rapidly exploding statistical 
error. We illustrate this with Fig. 2, where Monte Carlo simulations of the 
photon counting experiment for the one-photon Fock state are superimposed 
on the exact expectation value with the analytically calculated statistical error. 
Detailed discussion of these issues can be found in Ref. 3. 
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